Perturbation theory for optical excitations in the 
one-dimensional extended Peierls— Hubbard model 

Anja Grage, Florian Gebhard and Jorg Rissler 

Fachbereich Physik, Philipps-Universitat Marburg, D-35032 Marburg, Germany 

Abstract. For the one-dimensional, extended Pcicrls-Hubbard model we cal- 
culate analytically the ground-state energy and the single-particle gap to second 
order in the Coulomb interaction for a given lattice dimerization. The comparison 
with numerically exact data from the Density-Matrix Rcnormalization Group 
shows that the ground-state energy is quantitatively reliable for Coulomb 
parameters as large as the band width. The single-particle gap can almost 
triple from its bare Peierls value before substantial deviations appear. For the 
calculation of the dominant optical excitations, we follow two approaches. In 
Wannier theory, we perturb the Wannier exciton states to second order. In two- 
step perturbation theory, similar in spirit to the GW-BSE approach, we form 
cxcitons from dressed electron-hole excitations. We find the Wannier approach 
to be superior to the two-step perturbation theory. For singlet excitons, Wannier 
theory is applicable up to Coulomb parameters as large as half band width. For 
triplet excitons, second-order perturbation theory quickly fails completely. 



PACS numbers: 71.10Fd,71.35.-y, 71.27. +a 
1. Introduction 

Electrons in low-dimensional systems behave very differently from their three- 
dimensional counterparts. On the one hand, at low temperatures they tend to distort 
the lattice (Peierls effect For example, poly-acetylene as the simplest 7r-conjugated 
polymer shows an alternation between short and long bonds At the Peierls tran- 
sition a gap opens in the single-particle spectrum and the metal turns into a Peierls 
insulator. On the other hand, the electron-electron interaction provides a relevant 
perturbation in one dimension. For example, metals in one dimension are Luttinger 
liquids in which the low-energy charge and spin excitations propagate with different 
velocities so that quasi-particles as in Landau Fermi-Liquid Theory are absent in the 
single-particle spectral function [HIE]- At commensurate band fillings, the electron- 
electron interaction induces Umklapp scatterings which turn the metal state into a 
Mott insulator. 

This immediately raises the question: which of the two perturbations, Peierls 
distortion or Mott transition, dominates the other one? This question can be 
answered within the framework of field theory because both the electron-lattice and 
the electron-electron interactions destroy the metallic phase even in the limit of weak 
coupling. As has been shown in E|, the Peierls distortion is always a relevant 
perturbation in the sense of a renormalization-group analysis whereas the Umklapp 
scatterings which cause the Mott transition are only marginally relevant. At large 
Coulomb interactions when the system is a Mott insulator, the Peierls distortion is also 
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present at low temperatures because the effective spin system undergoes a spin-Pcicrls 
transition. Therefore, one-dimensional lattices are expected to be Peierls distorted at 
zero temperature, irrespective of the presence of the electron-electron interaction. 

In presence of a Peierls distortion and one electron per lattice site (half band- 
filling), the Peierls insulator can be described in terms of a filled lower Peierls band 
and an empty upper Peierls band which are separated by a finite Peierls gap whose 
size is proportional to the dimerization strength. In principle, the finite Peierls gap 
permits a perturbativc treatment of the electron-electron interaction. It should be 
clear, however, that such a perturbative treatment is limited to 'weak coupling', 
and the range of the validity of perturbation theory remains to be determined. 
In particular, even the screened electron-electron interaction in polymers is quite 
sizable |7j, and it remains to be studied in how far Wannier theory and similar 
perturbative approaches [HI Ej lead to a reliable description of optical excitations in 
these materials. 

In this work, we investigate the extended Peierls-Hubbard model at half band- 
filling which models the Peierls distortion and the Coulomb interaction in ideal poly- 
acetylene. The model is introduced in section [3 In section |2| we use standard Ray- 
lcigh-Schrodinger perturbation theory to calculate the ground-state energy to second 
order in the Coulomb interaction and compare our results to numerically exact data 
from the Density-Matrix Renormalization Group (DMRG). We show that second- 
order perturbation theory is applicable up to interaction strengths as large as the 
band width. In section 0] we observe the same behavior for the gap for single-particle 
excitations which increases quickly as a function of the interaction strength. 

In order to study optical excitations, we follow two routes because Rayleigh- 
Schrodinger perturbation theory is inapplicable for a description of bound electron- 
hole pairs. In scction[S]we develop Wannier perturbation theory to second order where 
an exciton is formed already in first order. Alternatively, we introduce the 'two-step 
perturbation theory' in which the exciton is formed in a second step after the electrons 
and holes have been dressed first. In sectionElthc comparison of our analytical results 
with those from the DMRG show that second-order Wannier theory is superior to 
the two-step perturbation theory. The singlet exciton is reasonably well described by 
Wannier theory for moderate interactions but perturbation theory quickly fails for 
the triplet exciton. We close our presentation with concluding remarks in section 
Technical details of the calculations are deferred to the appendix. 

2. Hamilton Operator 

2.1. Peierls Model 

We investigate spin-1/2 electrons on a periodically distorted chain of L sites whose 
motion is described by the Peierls Hamiltonian, 

f = (-t) £ (1 + (-1)'*) (c+ a c l+l a + c+ +l a Q , (1) 

where ct . ci )C are creation and annihilation operators for electrons with spin a =|, | 
on site I. The lattice spacing is set to unity, ao = 1. Since we are interested in 
the insulating phase, we consider exclusively a half-filled band where the number of 
electrons N = 2N a equals the number of lattice sites L. For our analytical calculations, 
we choose periodic boundary conditions and L/2 even. For the DMRG investigations, 
open boundary conditions are employed. 
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The Peierls operator is diagonal in momentum space |l(Jj . 

f = Y J E(k)(bUK.-^k,a) , (2) 

k-.cr 

where k = 2nm/L 1 m = — L/2, . . . , L/2 — 1, are the momenta from the reduced 
Brillouin zone. The dispersion relation for the upper and lower Peierls bands is 

E{k) = Ve(fc) 2 + A(/c) 2 (3) 

with 

e(k) = - 2icos(fc) , (4) 
A(k) = 2t8sm(k) . 

The Fermi operators for the electrons in the Peierls bands obey 

Sfe.cr = a>kCk,a + i/SfcCfc+Tr^cr , 

h,(T = hc-k,a - iafcC fc +7T,<T (5) 

with 



1 (a W\ 

2 { E(k)) ' 



and sgn(a; 7^ 0) = x/\x\ is the sign function. For the inverse transformation and other 
useful relations, see |Appendix A| 



2.2. Coulomb interaction 

The electrons arc supposed to interact locally with strength U. The Hubbard inter- 
action reads 

L> = C/^(n i . T -l/2)(n u -l/2), (7) 
l 

where hi !<7 = c/" CT c ; a is the local density operator at site I for spin a. Screening is not 
very efficient in an insulator. Therefore, we take into account the long-range Coulomb 
interaction of the electrons in the form of a \/r potential of effective strength V, 

V^VW'Efa-Wt+r-l) , (8) 

r I 

where hi = hi^ + hi t i counts the electrons on site I and V(r) = V/(2r) describes the 
Coulomb potential. For the analytical calculations, the specific form of V(r) is not 
important. 



2.3. Extended Peierls-Hubbard model 

Altogether we investigate the one-dimensional extended Peierls-Hubbard model, 

H = f + U + V = f + W . (9) 

The Hamiltonian is invariant under SU(2) spin transformations so that the total spin 
is a good quantum number. Including the charge-SU(2) the Hamiltonian is SO(4) 
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symmetric. Moreover, it exhibits particle- hole symmetry, i.e., H is invariant under 
the transformation 

PH: c+ a ^(-l)% a ; cv (-1)'^ • (10) 

The chemical potential /i = guarantees a half- filled band for all temperatures 

The Pcierls dimcrization (<5 7^ 0) and the Coulomb interaction (U, V 7^ 0) 
individually lead to an insulating ground state at half band-filling. A field-theoretical 
investigation shows El that the lattice distortion is always a relevant perturbation, 
i.e., it occurs for all values of the Coulomb interaction. Therefore, the Peierls 
insulator provides a valid starting point for a perturbation expansion in the Coulomb 
interaction. 



3. Ground-state energy to second order 

3.1. Analytical results 

For the calculation of the ground-state energy E we apply standard Rayleigh- 
Schrodingcr perturbation theory to second order, 



(11) 



The expectation value of H with the ground state of the Peierls insulator (filled lower 
Pcierls band) 



E 


= ^ 0) 










- E o > 




P (2) 
^0 


= E? 







|FS) J vacuum) 



(12) 



contributes 



E^°> = (FS|f|FS) = -2j2^Vcos(k) 2 + S 2 sin(fc) 5 



E% = (FS|[/|FS) = 



(13) 

EX = <FS|F|FS) =-L Yl 1V{r)[{Mr)f + (B s (r)) 2 } , (14) 



-L/2<r<L/2 
r|3=l 



where (r|2 = 1) denotes all odd r, ^(r), Bs(r) are defined in (| A. 5|) . and (|A.10|) is used 
for the derivation of (|13|) . For the Peierls insulator we have in the thermodynamic 
limit, 



lim 



(0) 



-(4t/7r)E(l-6 2 ) 



L — >OG L 

where E(m) is the complete elliptic integral of the second kind. 



(15) 



3.1.1. Second-order excitations. To second order, particle-hole excitations of the 
Peierls ground states must be considered. In general, we denote excitations with 
r (s) particle-hole excitations in the "["-sector (|-sector) by \rs). Fig. ^ shows two 
particle-hole excitations with antiparallel spins, i.e., |11), and parallel spins, i.e., |20), 
respectively. These are all excitations which contribute to the ground-state energy to 
second order. 




3.1.2. Second order in the Hubbard interaction. For the Hubbard interaction, only 
the states |11) contribute. As input to l(TT|l they give [5] 

|(FS|L>|11)' 2 



E 



p(o) _ p (o) 

|11) ^0 ^111) 



E 



ki ,&2 ,&4 



E(h) + E{k 2 ) + E(k 3 ) + E{k i ) 



|<5fei-fe2+fc3-fc4,o["l(fcl, fc 2, ^3, 
+ 4 1 -fc 2 +fc 3 -fc 4 ,±7r['"2(fcl, ^2,^3, ^4)] 2 } (16) 



with ^1,2(^1, fe, ^3, ^4) from ljA.4l 



3.1.3. Second order in the long-range interaction. The long-range Coulomb inter- 
action also induces scatterings between particles with the same spin. Therefore, we 
have three contributions to second order. We find jS] 

1 . 2 



K |iU/ =E 



E 



10) ^0 

where vo(k) is given in (|A.6|I 
,y 2 |ii) _ 



(FS|V|10) 



(0) 



E, 



(0) 

1 10) 



E 



V 2 v (k) 2 
E(k) 



E> 



E 



fci ,k2 



S(fci) + £(fc 2 ) + £(fc 3 ) + E(fol) 
|<Jfci-fca+*3-Jfe4,obl(Axi fe, &3, ^4)] 2 
+ 5fc 1 -fea+fc 3 -fe4,±7r[v2(fcl) ^2,^3, ^)] 2 } 



(17) 



(18) 
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with v\ t 2{k\, k 2 , &3, ^4) from (|A.8|I and (|A.9|I . and 

2 



V 2 |20) 



E 



E(h) + E{k 2 ) + E(k 3 ) + E{k±) 



5k 1 -k 2 +k 3 -k i ,a[vi{ki,k 2 , h, k 4 ) - vi(ki, k A , fc 3 , k 2 )} 2 

+ Ski-k 2 +k 3 -k A ,±n[v2(ki, k 2 , fc 3 , ki) - v 2 (k\,ki, fc 3 , k 2 )] 2 . 

(19) 

The summation restriction in the last term can be relaxed. Due to the symmetry 
of the functions ^1^(^1,^2,^:3,^4) we may multiply the last contribution by a factor 
1/4 and sum over all momenta from the reduced Brillouin zone. Due to the spin-flip 
symmetry we altogether have 



EX 



2E, 



v 2 |io) 




E 



V 2 |ll> 





2E, 



y 2 |2o) 




(20) 



as input to Qllfl. 



3.1.4- Second-order mixed interactions. Only the states |11) contribute to order UV, 



E, 



uv 



(FS|V r |ll)(ll|C/|FS) +c.c. 



p (0) p (0) 



E 



1 



E 

111) 

2UV 

k^tk.M E ^ + E ^ + E ^ + E ^ 
Sk 1 -k 2 +k 3 -k i ,oUi(ki,k 2 , fc 3 , fc 4 )ui(fci, k 2 ,k 3 , k A ) 

+ fik 1 -k 2 +k 3 -k 1 ,±TrU 2 (ki, k 2 , fc 3 , k 4 )v 2 (ki, k 2 ,k 3 , fc 4 ) 



(21) 



3.2. Comparison with numerical results 

3.2.1. Finite-size effects. First, we demonstrate that lattices with L = 100 provide 
results which are very close to the thermodynamic limit. Fig. gives an example for 
U = 2V = 2t and S = 0.2. 

The result from perturbation theory for periodic boundary conditions shows 
almost no size dependence whereas the data from the numerical density-matrix 
renormalization group displays the typical parabolic convergence in 1/L. The 
comparison shows that the results for L = 100 are almost identical to the result in the 
thermodynamic limit. In particular, the differences between perturbation theory and 
the numerically exact DMRG are small but significant. The same observation holds 
equally well for other choices of (finite) model parameters U, V, 8. 

In the following we shall show results for L = 100 and take them as representative 
for the thermodynamic limit. 



3.2.2. Fixed ratio U/V. Next, we discuss the ground-state energy density as a 
function of the Coulomb interaction for a fixed ratio of U /V. This is shown in 
Fig. |3| where we compare the DMRG data with the results from first and second- 
order perturbation theory It is seen that the second-order correction considerably 
improves the agreement between the numerically exact DMRG data and perturbation 
theory. 
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Figure 2. Ground-state energy density Eq/L as a function of inverse system size 
for fixed U = 2V = 2t and <5 = 0.2 in the DMRG and perturbation theory; lines 
are polynomial fits. Note the resolution of the ordinate. 
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Figure 3. Ground-state energy density Eq/L as a function of V for fixed 
U/V = 2, L = 100, and S = 0.2 in the DMRG (circles) and perturbation theory. 
Lines are polynomial fits. For U = V = 0, lim^—,^ Ej^' /L m 1.34. 
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3.2.3. Fixed V . In general, the influence of the Hubbard interaction is well described 
by the second-order term, as can be seen from Fig. It shows that perturbation 
theory works very well even when U is as large as the band width W = At. The 
agreement is even better for smaller V. 




u/v 

Figure 4. Ground-state energy density Eq/L as a function of U for V = 2t, 
L = 100, and <5 = 0.2 in the DMRG (circles) and perturbation theory. Lines are 
polynomial fits. 



3.2.4- Conclusion. The above comparison has been made for a sizable value of 
the Peierls dimerization, 5 = 0.2, so that the bare Peierls gap, A p = 4t5 = 0.8t, 
is substantial. Naturally, the comparison between the DMRG and second-order 
perturbation theory becomes less favorable when we decrease 6 to realistic values, 
5 < 0.1; see section [7| Nevertheless, even for small values of 5, the formulae in this 
section give a very good estimate for the ground-state energy density for a large range 
of Coulomb interactions, U « 2V < W. 

4. Single-particle gap to second order 

In order to generate current-carrying excitations in a Peierls-Mott insulator, the gap 
for single-particle excitations must be overcome. It is defined by 

Ai = [E (L + 1) - E (L)] - [E (L) - E (L - 1)] , (22) 

where Eo(N) is the ground-state energy of the system with TV electrons. Due to 
particle-hole symmetry we have 

Ai = 2[E (L + 1) - E (L)] . (23) 
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4-1- Analytical results 

For the perturbative calculation of the single-particle gap Ai we apply standard 
Rayleigh-Schrodinger perturbation theory again and write 



.(<» 



A^ + A{ +A^+ AJ 7 + A^ 



r/v 



(24) 



For the Peierls insulator the first two contributions are readily calculated because the 
ground state with N = L + 1 particles is 

|p)^6+ T |FS) , P=-tt/2, (25) 

and, thus, 



(o) _ 



(p\f\p) - (FS|T|FS) = 2E( P ) = US = A p 



(26) 



is the Peierls gap for all system sizes. Due to the spin symmetry we were allowed to 
use a =|- Moreover, the first-order contributions read 

, (27) 



Af 

Ar 



(p\U\ P ) - (FS|f7|FS) 
2 [<p|V» - (FS|t>|FS) 

2 4 * y ( r ) 



cos(p) cos(pr)As(r) + 5sin(p) sin(pr)Bg(r) 



-L/2<t<L/2 
r|2=l 



E(p) 



Y, 4^(r)(-l)[ (r - 1 )/ 2 ] J B 5 (r 



(28) 



-L/2<r<L/2 
r|2=l 



where p = —tt/2 was used in the last step. 



4.1.1. Second order in the Hubbard interaction. In the intermediate states an 
additional j electron and r (s) electron-holes excitations in the t-sector (|-sector) 
are possible, denoted as \k;rs). There are two contributions to second order in the 
Hubbard interaction, 

Af |01> =2 ]T 



and 



|(p|CT+ a a fa)i |fc 3 )| i 



ki,k 2 ,k 3 
2 



E(p) - (E(h) + E(k 2 ) + E(k 3 )) 

2 (1 - Spjej) 



E 



Ai 



U 2 \ll) 



E 



E(p) - (E(ki) + E(k 2 ) + E(k 3 )) 
[4i-fa+fa-p,oh(^^2, h,p)] 2 

+ h 1 -k3+k 3 -p,±A u lO ! U fa, k 3 ,p)} 2 

\(j>\Ub+ iA a k2!l b+ 3il a kiil \p)\ 2 



2 



E 



(E(h) + E(k 2 ) + E(k 3 ) + E(k 4 )) E ° 
5 P ,ki 



(29) 



t/ 2 |ll> 



k± .k% 7 k% ,/c4 



E{k l ) + E{k 2 ) + E(k 3 ) + E{k A ) 
^ fel -fc 2 +fc 3 -fc4,o[ u i( fc i> fc 2, k 3 , fc 4 )] 2 (30) 
+ fik 1 -k 2 +k 3 ~k i ,±A u 2(ki, k 2 , k 3 , ki) 
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(31) 



as input to l)24[l. 



4-1-2. Second order in the long-range interaction. To second order in the long-range 
interaction, the states |fc; 10) and \k; 01) contribute. Note that the fermionic nature of 
the operators must be taken into account properly, e.g., b^\p) = — ||fc). This gives 



the contributions 



A 



y 2 |io) 



E 



2|<p|vS+ iT a fc2 , T \k 3 )\< 



k 2 ,k 1 >k 3 
2 



v t.) E 



E(p) - (E(h) + E{k 2 ) + E{k 3 )) 
2 



2E, 



v 2 \w) 
o 



k 2 ,k\>k 3 



E(p) - (E(h) + E(k 2 ) + E(k 3 )) 



0ki — k 2 + k 3 —p,0 ^S p ^k 3 Sk u k 2 Lv o(ki) - & p , kl 8 k2 ^ 3 LvQ(k 3 ) 

+ v 2 (k 3l p, k 1} k 2 ) - v 2 (k 1: p, k 3 , fc 2 ) j 
+ <5fc 1 -fc 2 +fc 3 - P ,± 7r {vi(fc3,.P, h, k 2 ) - vt(ki,p, k 3 , fc 2 ) j 



2E, 



v 2 \io) 





(32) 



and 



A 



y 2 |oi) 



E 

k\ ,k2 ,&3 
2 



\(p\Vbl A a k2d \k 3 )\< 



E(p) - (E(h) + E(k 2 ) + E(k 3 )) Ea 
1 



y 2 |oi) 



k\ ,k 2 ,k 3 



E(p) - (E(h) + E(k 2 ) + E(k 3 )) 



Sk 1 -k 2 +k3-p,o\s P! k 3 6k 1 ,k 2 Lv (ki) + v 2 (k 3 ,p, ki, fc 2 ) j 

— k 2 +k 3 —p,±ir |«l(fc 3 ,p, fci,fc 2 )| 



V 2 |01) 



(33) 



Moreover, the states with two particle- hole excitations, \p; 11), contribute. We find 



a; 



v 2 ih> 



= E - 



\(p\Vb+ i1 a k2 /b+ al a ki . l \p)\' 2 



kl ,/C2 ,/C3 ,&4 

2 



E 



(E(h) + E(k 2 ) + E(k 3 ) + E(k 4 )) E ° 



V 2 \ll) 



k\,k 2 ,k 3 ,k& 



E(ki) + E{k 2 ) + E{k 3 ) + E(ki) 

[4i-fa+fa-p,oh(fci, k 2 ,k 3 ,ki)} 2 (34) 
+ S kl ^k 2 +k 3 - P ,±A v 2(ki, k 2 , k 3 , fc 4 )] J 
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and 



a; 



V' 2 |20) 



E 



1(^1^+ T a fc2 , T S+ ;T a fc4 , T |p)| 2 
^ - {E{k x ) + E(k 2 ) + E(k 3 ) + E{k 4 )) 



V 2 |20) 



E 



S(fc 1 ) + S(fe) + ^(fc 3 ) + J B(fc4) 



4i-fe 2 +A;3-p,o [vi(/ti, k 2 , k 3l k 4 ) - Vi(ki,ki, k 3 ,k 2 )] 2 
+ S kl -k2+k 3 -p,±A v 2{ki,k2, k 3l k 4 ) - v 2 (ki, ki, k 3 , k 2 )f 



with A j 



V 2 \02) 



Ap 20> . Altogether, from J^, §3, and we find 



l 10 > _L l 01 > _L A V 2 |ll> 



2A 



y 2 |2o) 



(35) 



(36) 



as input to (|24|l . 



4-. 1.3. Second-order mixed interactions. As for the Hubbard interaction we have two 
contributions for the UV mixed interactions. From the intermediate states |fc;01) we 
obtain 

2 



A 



r/r|oi> 



[<p|^£a a *».il*3)(p|v , &fc 1 ,i 6 fa.il*3>* + c - c - 



(^) 



j ^2 ) ^3 



S(p)-(£?(fei)+B(fe2) + S(fcs)) 



Ski-k 2 +k 3 -p,0V2(k 3 ,P, kl, k 2 )u 2 (kl,k 2 ,k 3 ,p) 

+ 5jn-fe 2 +fe 3 -P,±7rWi(fc3,p ! fci, k 2 )ui(kx,k 2 ,k 3 ,p) 



(37) 



and the states |p; 11) give 
Af y|11) = 



E 



(-2) 



E(h) + E{k 2 ) + E(k 3 ) + E(k 4 ) 



+ c.c] - E\ 



I7V]11) 



S(fc 1 )+^(fc 2 )+S(fc 3 )+S(fc4) 



fel 

4i-fc 2 +fc 3 -p,0Wl(fcl, fe, ^3, ^4)^1(^1,^2, &3, ki) 

+ 4 1 -fe 2 +fc 3 -p,±7rM2(fci, fe, £3, k 4 )v 2 (ki, k 2 , k 3 , ki) . (38) 

Altogether we thus find from these two equations 

A uv = A uv\oi) + A uv\ii) (3Q) 

as input to (J2J. 
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4-2. Comparison with numerical results 

4-2.1. Finite-size effects. As for the ground-state energy we first investigate the 
size dependence. Fig. [S] shows that the analytical results for the single-particle gap 
are almost independent of the system size, as expected. The DMRG data show the 
typical quadratic convergence as a function of 1/L. For our assessment of the quality 
of the perturbative approach we note that the DMRG data for L = 100 represent the 
value in the thermodynamic limit to an accuracy of about one percent. Therefore, the 
differences between perturbation theory and the numerically exact DMRG which we 
will discuss next are significant. 
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Figure 5. Single-particle gap Ai as a function of the inverse system size for 
fixed U = 2V = 2t and 5 = 0.2 in the DMRG and perturbation theory; lines are 
polynomial fits. 



4.2.2. Fixed ratio U /V . Fig. © compares the single-particle gap as a function of V 
for fixed ratio U/V = 2 in the DMRG and perturbation theory to first and second 
order. The agreement is worse than for the ground-state energy. Nevertheless, the 
deviation between the numerically exact DMRG results and perturbation theory 
at U = 3t and V = 1.5i is only about ten percent despite the fact that the 
gap has increased from its noninteracting value A p = 0.8t by a factor of 2.5 to 
A° MRG ([/ = 3t, V = 1.5t, S = 0.2) » 2M. 

The figure also shows that the second-order contribution only slightly improves 
the first-order result. The surprisingly good first-order result is partially due to the 
error compensation in the calculation of the ground-state energies for N — L and 
N = L + 1 particles. 
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V=U/2 



Figure 6. Single-particle gap Aj as a function of V for fixed U/V = 2, L = 100, 
and <5 = 0.2 in the DMRG (circles) and perturbation theory. Lower curve: 
perturbation theory to first order (PT1), upper curve: perturbation theory to 
second order (PT). Lines are polynomial fits. 




U/V 



Figure 7. Single-particle gap Ai as a function of U/V for V = t, L = 100, and 
5 = 0.2 in the DMRG (circles) and perturbation theory; lines are polynomial fits. 
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4-2.3. Fixed V . As in the case of the ground-state energy, the discrepancies between 
the DMRG single-particle gap and the result from second-order perturbation theory 
are mainly caused by the long-range part of the Coulomb interaction. This is shown 
in Fig. [7| The results from second-order perturbation theory closely follow the DMRG 
data for all U/t at fixed V = t. 

At = 3* and V = t the difference between A? MRG (f7 = 3*, V = t, 6 = 0.2) = 
1.8680* and A^ T (U = 3t,V = t,S = 0.2) = 1.7229* is only eight percent. Here, first- 
order perturbation theory gives A^ T1 (U = 3f, V = t, 5 = 0.2) = 1.430* which is off 
by 23 percent. This shows that the inclusion of the second-order terms improves the 
agreement considerably for all U/t as long as V/t is not too large. 



4-2.4- Conclusion. Second-order perturbation theory provides a reliable estimate for 
the single-particle gap in the regime where it has almost tripled from its bare Peierls 
value, i.e., in the region (U < At, V < 2t) for S = 0.2* (A p = 0.8*). This is rather 
important because it indicates that only a minor part of the single-particle gap in 
polymers is caused by the Peierls effect. For poly-acetylene where * f=a 2eV and 
Ai « 1.6 eV one would obtain S = 0.2 in the absence of Coulomb correlations. When 
we assume a correlation enhancement for the gap by a factor of about three, we would 
estimate S « 0.07, in agreement with other estimates in the literature. 

A reliable estimate for the single-particle gap is important for the calculation of 
binding energies of bound particle- hole excitations (excitons). They are the subject 
of the next section. 



5. Perturbative approaches to optical excitations 

There are various possibilities to set up second-order perturbation theory for optical 
excitations. In this work we discuss Wannier, two-step, and down-folding perturbation 
theories separately before we calculate matrix elements. Starting point is the subspace 
of optical excitations in the Peierls model, 

singlet : \s) = ^ w s (k)\k) s , 

k 

triplet : \t) = ^w t {k)\k) t , (40) 

fc 

with the normalization ^2k[w( s ,t)(k)] 2 = 1 and 

= Vf ()fc,T«M ± bt,i^,l) IPS) . (41) 

Here, we introduce linear combinations of single particle-hole excitations 1 10} and 1 01) 
from the outset because we are mainly interested in bound states of particle and hole. 



5.1. Wannier perturbation theory 



Let P and Q be projectors onto separate subspaces of the Hilbcrt space, P + Q = 1. 
In the case of excitons we let P project onto the subspace of a single particle-hole 
excitation, P = P±. Then, we may write for H = T + W 



H 



T + PWP 



QWP + PWQ + QWQ 



= H + H ± , (42) 
where Hq describes the unperturbed system and H± represents the perturbation. 



Optical excitations in the one- dimensional extended Peierls- Hubbard model 



15 



5.1.1. Wannier perturbation theory to first order. In the usual Wannicr theory, one 
diagonalizcs H in the subspacc of a single particle-hole excitation, 

H \(s,t)) n = et t) \(s,t)) n , (43) 

where e„ ' ; is the Wannier spectrum in the spin singlet and triplet sector. Typically, 
the states with lowest energy correspond to bound states (excitons) |12j . 

5.1.2. Wannier perturbation theory to second order. Given the Wannier spectrum 
and states to first order, we can calculate the Coulomb corrections to the Wannier 
eigenstates via perturbation theory to second order in H± 

rp \ - (x n \H±\m)(m\H ± \x n ) 

E n =(-n + }^ — (0) ( 44 ) 

\rn) £ ™ ^ m 

in the singlet and triplet sector, x = (s,t). Here, E$ = (m\Ho\m), and \m) can be 
chosen as an eigenstate of the Peierls Hamiltonian T because corrections are of higher 
order in the Coulomb interaction. The states |m) contain up to three particle-hole 
excitations. The corresponding matrix elements in (|44f) are calculated in the next 
section. 



5.2. Two-step perturbation theory 

The Wannier theory for excitons probes how a bound excitation of bare particles 
and holes is perturbed by the Coulomb interaction. In two-step perturbation theory 
the two steps are reversed: excitons are formed from dressed particles and holes. 
This two-step approach is very similar in spirit to the GW-BSE approach where the 
Bcthe-Salpeter Equation is solved for quasi-particles which are calculated in the GW- 
approximation j^j. 

To set up the two-step perturbation theory in general, we split the Coulomb 
interaction into two parts, W = Va + Vb- Let P r be the projection operators onto 
subspaces with r particle-hole pairs and Q r = 1 — P r its complement. We identify 

V A = PrWQr + QrWPr + Q r WQ r , (45) 
V B = PrWPr ■ (46) 

In the first step of the perturbation theory we calculate the change of an eigenstate \n°) 
of T due to the influence of the perturbation Va ■ In Raylcigh-Schrodinger perturbation 
theory for Ha — T +Va to first order we find 

l„ ) = z L. ) + E MMVdi<>. (47) 

( \m) ~ ^ m 

To second order the normalization factor Z n reads 

\n°\V A \m) 2 



\Z n \ 2 = l-^ 



(48) 



The states \n) are so-called 'dressed states' because they contain particle-hole 
excitations which are renormalized by the interaction Va- 
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In the second step we apply (almost degenerate) perturbation theory for the 
perturbation Vb- Therefore, we diagonalizc H = Ha + Vb in the basis of the states 
|n), i.e., we calculate the eigenvalues of the Hamilton matrix 

H{ ni ,n 2 ) = {ni\H\n 2 ) (49) 

in the respective subspaces with r particle-hole excitations. 

As an example we consider the case r = 0, i.e., we calculate the second-order 
ground-state energy in two-step perturbation theory. Since |0) is not degenerate, we 
have 

10°) = |FS) , (50) 
|0) = Z (\FS) + \X}) . (51) 
The Hamilton matrix is a scalar, 
E = (0\H\0) 

= (\Z \ 2 - 1) (FS|T|FS) + (FS|f + H^|FS) 
+ (X\W\FS) + (FS\W\X) + (X\f\X) 

\{FS\W\m) 2 

= (FS\f + W\FS) + J2 '(„)„«,) • ( 52 ) 

| m )^|FS) ^0 ^to 

This expression is identical to Rayleigh-Schrodinger perturbation theory to second 
order. 

5.3. Down-folding perturbation theory 

5.3.1. Equivalence to Brillouin-Wigner perturbation theory. The so-called 'down- 
folding' approach is also based on a projection to relevant degrees of freedom. Consider 
the stationary Schrodinger equation at energy E 

H\V) = (T + W) I*) =.E7|*). (53) 

For any two projectors P and Q with P + Q = 1 we define the orthogonal states 

\<p)=P\9) , \X) = Q\*) (54) 
which obey the coupled Schrodinger equations 
E\tp) = PH\ V ) + PH\ X ) , 

E\ X ) = QH\<p) + QH\ X ) ■ (55) 
The second equation is formally inverted with the help of the resolvent operator 
(z — H)' 1 to give a single equation 

PHP\tp) + PHQ (e-H\ 1 QHP\ip) = E\tp) . (56) 

We now set P = P r so that 

P r HP r =f + P r WP r , P r HQ r = P r WQ r , Q r HP r = Q r WP r , (57) 
and insert the unit operator in (|56|l to obtain 

ry zj ry \ \ i PrWQ r \m) (m\Q r W r P r \^ 

PrHP r \cp) +2_^ F _ F =E lp \<p). (58) 

Hits) I-Jm 

\m) 

The same expression is obtained in Brillouin Wigncr perturbation theory in which 
the energy E = E v must be determined self-consistently. 



Optical excitations in the one- dimensional extended Peierls-Hubbard model 17 

5.3.2. Size-consistent reformulation. The down-folding or Brillouin Wigner pertur- 
bation expansion is hampered by the fact that it is not size-consistent. Therefore, the 
form Q58f) cannot be used for our calculations. 

In order to make progress, we first subtract the ground-state energy from all 
energies in (l58l) . 

(t? bmi \ t>(tr w\\ y , \- P r WQ r \rn)(m\Q r W\<p) 

so that all energy differences correspond to excitation energies of order unity. To 
second order in the interaction we may write 

E v - E = e v , E m -E = E^ - + 0(U, V) . (60) 
Therefore, to second order we find 

I \ Pffr ni y , \- P r WQ r \m)(m\Q r W\y) 
e» =P r (H-E Q M+^ -— - M — . (61) 

\ m ) e v> ^0 ) 

As an example, we apply this equation for the ground state, i.e., r = 0, \ip) = |FS), 
and e v = 0. When we multiply 1|61|) from the left with (FS| we immediately recover 
the result from Rayleigh-Schrodinger perturbation theory. 

For the exciton states \x) and energies e x we find the Schrodinger equation 

l\ Pm TP m v i V- PiWQi\m)(m\QiW\x) 
e x \x)=P 1 (H-E )\x)+}^ — Z(o)— • (62) 

\m) ex ~ ~ ^0 I 

This translates into a matrix-diagonalization problem with the diagonal entries 

H BW (e (Stt y,k u ki) = is , t) (ki\f\ki) M - (FS|f |FS) - E™ 

+ " (FS|#|FS) 

, ( s ,t) I T^|m> (s , t) (Sjt) <m| VT| fei > (S;t) 
+ E " A F (o), (6 3 ) 

\m) e (M) _ K^rn - E Q ) 

and the non-diagonal entries (ki 7^ k%) 

HBw(e( s ,t);k 1 ,k 2 ) = (s,t){ki\W\k 2 )( s j) 

(s,t){ki\W\m)( S}t )(^ t) {m\W\k 2 )( s .t) 
|m) e (M) ~ ^ m ~ h o > 

(64) 

where the excited states |m) contain none, two, or three particle-hole excitations of 
the Peierls ground state. Note that the eigenvalue e( S)t ) of the matrix i?Bw(e( s ,t)) 
in (|64() must be determined self-consistently. In our calculations we target the lowest 
eigenvalue as the bound exciton in the triplet and singlet sectors. 

6. Excitons to second order 

6.1. Analytical results 

The calculation of optical excitation energies requires the matrix elements 

M ^ Stt )(pi,P2) = (s,t)(Pl\T\P2)(s,t) ~ S P1,P2 E 0) = S Pl,P2 2E (Pl) ! ( 65 ) 
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and 

M 1:(s ^(pi,p 2 ) = { s ,t){Pi\W\p 2 )(sd) ~ Spi, P 2 E o > ( 66 ) 
(s,t)(jPi\ W\m) (s j )[s , t) {Tn\W \p 2 ){ s .t) 

where |m) contains up to three particle-hole excitations. 



w / x \- («,t) ffi I ^ \ m )(s,t) (m) Pa^.t) , . 



6.1.1. Matrix elements to first order. The matrix elements M\ are readily calculated. 
We find 



with 



and 



M 1)(M) (pi,p 2 ) = M^ (st) ( Pl ,p 2 ) +Ml (s !t) (pi,p 2 ) (68) 

^ (s , t )(Pl:f2)=±^, (69) 



„v , s T r/ \ r o cos \P\ j cos \ P\ r ) 4 . , 
M?, {s!t) (pi,P2) = 5 PuP2 }2 E( Pl ) As{r) 



-i/2<r<i/2 
r!2=l 



>cos(pi) cos(pir) 
%) 

8<5 sin(pi) sin(pir) 



^ ^V(r) cos((pi - p 2 )r) 

[/lbi,p 2 ) + (-i) r / 1 2 bi,P2)] 



+ 7E^)t 1 + ( 1±1 )(- 1 ) r ] • ( 7 °) 

' r=iO 



6.1.2. Matrix elements to second order. The intermediate states 1 00) = |FS), |11), 
|20 ) |02), 1 30), |0 3), 1 21), and [12) contribute to fgj). The respective terms are lengthy, 
see .Appendix B| As in sections [3] and 01 we write 

M 2 ,( s ,t)(e;pi,P2) = M^( st) (e;p 1 ,p 2 ) + M^' {st} (e;p ll p 2 ) 

+ M^ t) (e; Pl ,p 2 ) . (71) 

The contributions to M^ 8 t )(e;pi,p 2 ) can be found in ijB.lfl and (|B.2fl . The terms for 
M% 3 , a t) (e;pi,p2) result from l|B~3)l . (jBTj), (jR5j), (|B~7|) . and l|B~6)l . Lastly, the terms 
for McfK t )(e;pi,p 2 ) arc collected in (|B.8ll and (|B.9|) . 



6.1.3. Wannier perturbation theory. The first-order result from Wannicr pcrturba- 

-(s,t) 



tion theory for the singlet /triplet excitons follows from the lowest eigenvalue e\l\\ of 



the Wannier matrix with entries 

M w { Pl , P2 ) = M 0!(s , t) Oi,p 2 ) + M 1)(M) (pi,p 2 ) . (72) 

Moreover, this calculation gives the first-order Wannier wave functions with (real) 
amplitudes W( s . t )(p), cf. (f^Ufl . 
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To second order, we evaluate Q44JI. 

(2) ^ ((s,t)\H±\m){m\H ± \(s,t)) 

On the right-hand-side of this equation we keep only the kinetic energy of the Wannicr 
exciton. 

e S)=E 2S ( fc )KM)(fc)] 2 (74) 

k 

in order to be consistent to second order in the Coulomb interaction. The second-order 
correction of the Wannicr excitons' excitation energy is then given by 

e fs,t) = ]C w (s,t){Pi)w(si){P2)M 2 . ( ^ t) {e [ ^ t yp l , P2 ) - Ef> , (75) 

Pl,P2 

and the excitons' excitation energy is e.{ s> t) = e (s t) + e (s *)' 

6.1.4- Two-step perturbation theory. To first order, the results from the Wannier and 
two-step perturbation theories are identical. For the calculation of the corrections we 
define the matrix 



B(„,t)(puP2) = ^(s,t)(Pi\W\m)(m\W\p 2 )(s, 







P (0) p(0) 
Him — -C/q 



(2E( Pl ) - (E { ° ] - E^ 0) ))(2E(p 2 ) - (E^ - E^)) 



(76) 

In order to calculate the entries of this matrix the energy denominators which appear 
in | Appendix B| must be modified appropriately. 

As shown in Q49[) . the excitation energy of the excitons in two-step perturbation 
theory is then obtained from the lowest eigenvalues of the matrix with the entries 

A( s .t){pi,P2) = M 0> ( s j- ) (p 1 ,p 2 ) + Af li(s!t) (pi,p 2 ) 

+ S PltP2 (Af 2 , (att) (2B(pi);pi J pi) - 4 2) ) 

+ (1 - 5p I ,p 3 )(M2,( s ,t)(2J5(pi);pi ) p 2 ) 
+ M Us j ) (2E(p 2 ); Pl ,p 2 ) 

+ B {a , t) (pi,P2)) ■ (77) 

6.1.5. Down-folding perturbation theory. According to l|63|) and (|64|l the down- 
folding method requires the diagonalization of the matrix with the entries 

F(s,t)(e(s,ty,Pi,P2) = M 0> ( s rf(p 1 ,p 2 ) + M lt{s j ) (p 1 ,p 2 ) 

( 2") 

+ M 2j ( a ,t)(e( ait );pi,p2) - $pi. P 2 E o • (78) 
The lowest eigenvalue must be determined self-consistently. 
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6.1.6. Comparison of numerical effort. The numerically cheapest method is the two- 
step perturbation theory. For given Peierls dimcrization S and functional form of V(r) 
one can calculate and store the matrix elements which define ^4( s ,t)(pi,P2)- The 
parameters U, V only enter as multiplicative factors. 

This docs not hold for the Wannier approach anymore because the kinetic energy 
depends on the wave function of the Wannier cxciton in first order of the interaction. 
For a first guess we could replace e|°^ in (|74|l by 2E{—-k/2) because p = — 7r/2 is the 
dominant contribution in the excitonic wave functions. This speeds up the analysis as 
a function of the Coulomb parameters but it also reduces the quality of the results. 

The down-folding approach is the most costly of the three approaches in terms 
of computer-time. For given interaction parameters, the effort of the second-order 
Wannier theory has to be repeated some six to eight times until convergence is reached. 

6.2. Comparison with numerical results for the singlet exciton 

We start with a comparison of our analytical results for the singlet cxciton with those 
of the DMRG. For the latter we assume that the singlet cxciton is identical to the 
lowest- lying charge excitation at fixed particle number N = L. This assumption is 
justified in the presence of a singlet exciton. 

We do not present results from the down-folding approach. The self-consistency 
algorithm is stable but the results from the down-folding approach underestimate 
the excitonic excitation energies drastically. We do not consider the down-folding 
perturbation theory any further because it is numerically the most expensive and 
quantitatively the least successful of our perturbative approaches. 

6.2.1. Finite-size effects. As in previous sections we first investigate the size- 
dependence of the excitation energy. In Fig.|H|we compare the results from the DMRG 
with those from the various perturbative approaches which we discussed in section [5] 
The DMRG data show a very weak size dependence which is a clear signal for a bound 
state. The data PT and PT1 do not reproduce this 1/L dependence. In fact, they 
describe an unbound particle-hole excitation. The starting point for PT and PT1 is a 
single particle-hole excitation at momentum p = — 7r/2, and corrections are calculated 
in first and second-order Raylcigh-Schrodingcr perturbation theory, respectively. The 
results show that this approach is not applicable for the description of excitonic optical 
excitations, and we shall not pursue this approach any further. 

Wannier theory to first order (WPT1), to second order (WPT2), and the two- 
step perturbation theory show a very weak size dependence, characteristic for a bound 
state. As seen from the figure, the two-step perturbation theory does not improve the 
result from first-order Wannier perturbation theory, and the second-order Wannier 
theory is the best approximation to the DMRG data. 

In all cases, the results for L = 100 are a very good estimate for the results in the 
thermodynamic limit. The differences between the various perturbative approaches 
and the numerically exact DMRG are significant. 

6.2.2. Fixed ratio U /V. In Fig.EDwe show the results for the singlet excitation energy 
as a function of V for L = 100, S = 0.2, and fixed ratio U/V = 2. It is seen that 
second-order Wannier theory is naturally superior to the first order approximation 
and provides a rather good description. At (V = l.bt, U = 3t), when the charge 
gap has almost doubled from its non-interacting value A p , the difference between the 
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Figure 8. Singlet exciton energy (charge gap) A c as a function of inverse system 
size for fixed U = 2V = 2t and 8 = 0.2 in the DMRG and perturbation theory; 
see text for further details. 




V=U/2 



Figure 9. Singlet exciton energy (charge gap) A c as a function of V for fixed 
ratio U/V = 2, L = 100, and S = 0.2 in the DMRG and in various perturbativc 
approaches. 
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DMRG and second-order Wannier theory is less than ten percent. In contrast, the 
first-order Wannier theory is already off by 30 percent, and two-step perturbation 
theory becomes even worse than that. 

6.2.3. Fixed V . As for the ground-state energy and the single-particle gap. the 
discrepancies between the DMRG and second-order perturbation theory increase as 
a function of V but not so much as a function of U. In Fig. ^] we show the singlet 
excitation energy as a function of U for fixed V = t, L = 100, and 5 = 0.2. The two- 
step perturbation theory and second-order Wannier theory closely follow the almost 
quadratical increase of the charge gap as a function of U. The singlet excitation energy 
can triple from its bare value A p = Q.8t at (U = 5i, V = t), and still the deviations 
are only about ten percent. In contrast, first-order Wannier theory begins to deviate 
noticeably already at (U = 3t, V = t). 




Figure 10. Singlet exciton energy (charge gap) A c as a function of U/V for 
fixed V = t, L = 100, and S = 0.2 in the DMRG and in various perturbativc 
approaches. 



6.2.4- Conclusion. Wannier theory to second order consistently improves the results 
from first-order Wannier theory. It provides a quantitative description even for large 
Coulomb parameters when the charge gap has more than doubled from its bare Peierls 
value. 

The two-step perturbation theory can handle the influence of large Hubbard 
interactions but it fails quickly when the long-range parts of the Coulomb interactions 
become relevant. This is somewhat unfortunate because the two-step perturbation 
theory is numerically the cheapest of all second-order methods. 

The good performance of the second-order Wannier theory indicates that this 
approach describes the relevant objects appropriately. To first order in the Coulomb 
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interaction Wannier theory forms the charge neutral exciton from a negative (electron) 
and a positive charge excitation (hole). This object has the appropriate quantum 
numbers also for strong Coulomb interactions. Therefore, the residual Coulomb 
interactions, i.e, the coupling to the ground state and other particle-hole excitations, 
apparently perturb the exciton properly. In contrast, two-step perturbation theory 
builds excitons from dressed electrons and holes and thereby overestimates the 
repulsion against other states with electron-hole pairs which results in a too small 
energetic splitting between the exciton in two-step perturbation theory and the ground 
state. 

6. 3. Comparison with numerical results for the triplet exciton 

We continue with a comparison of our analytical results for the triplet exciton with 
those of the DMRG. For the latter we assume that the triplet exciton is identical to 
the lowest-lying spin excitation at fixed particle number N = L. This assumption 
is justified in the presence of a triplet exciton. Again, we do not present results 
from the down-folding approach because it fails quantitatively already for fairly small 
interaction strengths. 

6.3.1. Fixed ratio U/V. For the triplet exciton, the results from perturbation theory 
are much less reliable than for the singlet exciton. As can be seen from Fig. ^2 a 
quantitatively correct answer is provided by second-order Wannier theory only up to 
U = 2V = 2t, whereas first-order Wannier theory and two-step perturbation theory 
fail for even smaller Coulomb parameters. 



< 




V=U/2 

Figure 11. Triplet exciton energy (spin gap) A s as a function of V for fixed 
ratio U/V = 2, L = 100, and 8 = 0.2 in the DMRG and in various perturbative 
approaches. 
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6.3.2. Fixed V . Perturbation theory fails even qualitatively for moderate to large 
values of the Hubbard interaction. The DMRG exciton energy increases, flattens out, 
and finally turns into a 1/U behavior, as seen in Figs.llllandll2l This shows that the 
nature of the excitonic excitation changes drastically as a function of U/t: for large 
Coulomb interactions the triplet exciton is a bound state of a pair of two spins on 
singly occupied sites, i.e., it looses its character of bound charge excitations which is 
assumed in perturbation theory 



T 




Figure 12. Triplet exciton energy (spin gap) A s as a function of U/V for fixed 
V = t, L = 100, and <5 = 0.2 in the DMRG and in various perturbative approaches. 



6.3.3. Conclusion. Triplet excitons cannot be described reliably within low-order 
perturbation theory because the nature of the lowest triplet excitation quickly changes 
as a function of the interaction strength. For not too large Hubbard interactions, the 
triplet exciton looses its charge contents and is better viewed as a bound state of 
two spin- 1/2 excitations (spinons). The description of such a state lies beyond the 
possibilities of finite-order perturbation theory. 

7. Concluding remarks 

7.1. Summary 

In this work we performed second-order perturbation theory for the extended Peierls- 
Hubbard model for the ground-state energy and the single-particle gap. The agreement 
with the numerically exact DMRG data was very good up to Coulomb parameters 
where the bare Peierls gap has almost tripled. 

For bound optical excitations, the situation is less favorable. The singlet exciton 
energy is reasonably well described by second-order Wannier perturbation theory for 
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Coulomb parameters which lead to a doubling of the single-particle gap. The triplet 
exciton can be described perturbatively for weak coupling only. 

The second-order Wannicr theory proves to be superior to other second-order 
approaches. Its success may be related to the fact that after its formation the exciton 
is a charge neutral object to which other parts of the Hilbert space are coupled less 
intensely by the Coulomb interaction. The computationally most costly down-folding 
method performs worst and can be safely discarded. The two-step perturbation theory 
where excitons are formed from dressed particles is computationally the cheapest 
method. Unfortunately, it does not improve the first-order Wannier theory for the 
singlet exciton systematically. This indicates that the exciton should not be seen as a 
bound state of dressed quasi-particles. In order to investigate this question further, the 
GW-BSE approach |U ought to be applied directly to the extended Peierls-Hubbard 
model so that its predictions can be tested against the numerically exact DMRG 
results. 

7.2. Discussion 

In order to compare our results with those of previous work on poly-(di)acetylene, 
we note that we may write V = e 2 /(ea ) where a ~ 1.4A is the average distance 
between the carbon atoms, e is the electric charge and e is the dielectric constant of the 
surrounding medium. In jS] its value is given by EAbc ~ 5 which is much larger than the 
typical values for plastic, e w 2.3. For the Abe value one finds V ~ 2 eV = t. and the 
application of first-order (Wannier) perturbation theory at 8 — 0.2 gives a reasonable 
agreement with the measured binding energies of the first singlet and triplet exciton 
in poly-diacetylene which are A c = 0.5 eV and A s = 0.9 eV, respectively. However, 
the measured single-particle gap is Ax = 2.3 eV which is seriously overestimated by 
0.7 eV for the given choice of parameters, Af MRG {U = 2t,V = t,S = 0.2) = 3.0 eV. 

The reason for this discrepancy is readily identified. The assumption 6 = 0.2 leads 
to a bare gap of Ap = 0.8i w 1.6 eV. As has been seen in Fig. El the single-particle gap 
increases linearly with the Coulomb interaction. Therefore, a good agreement with 
the experimental gap value is achieved at V = 0.6t which would correspond to e w 8 
which is unrealistic. Apparently, the assumption of a large Peierls gap is not correct. 

A more realistic approach starts from a much smaller Peierls contribution to the 
gap. We should rather work with S = 0.1 or smaller ^2], which is in accordance 
with an estimate of the electron-lattice interaction in benzene jjj. The Peierls 
gap then is Ap = OAt f» 0.8 eV which is about one third of the full gap. We 
then find for the single-particle gap A PT ([/ = 2t, V = t,S = 0.1) = 1.75 eV 
from second-order perturbation theory which is 12 percent below the exact value 
of A 1 3MRG (C7 = 2t,V = t,S = 0.1) = 1.98 eV from the DMRG. For the singlet exciton, 
we find A PT (L/ = 2t, V = t, 5 = 0.1) = 1.44 eV which is 22 percent above the DMRG 
value Aj? MRG (t/ = 2t, V = t, S = 0.1) = 1.18 eV. The agreement is quite acceptable in 
both cases. However, the comparison of the binding energies for the singlet exciton, 
A b£d, s ( U = 2t,V = t,6 = 0.1) = 0.31 eV versus A™f d RG ([/ = 2t,V = t,S = 0.1) = 
0.80 eV reveals that we are stretching second-order perturbation theory beyond its 
limits because the error in the binding energy increases to more than 60 percent. 

A realistic description of the experimental data for poly-diacetylene requires the 
inclusion of the full band-structure and of the lattice relaxation for excited states, 
mostly for the triplet exciton |14|. In particular, a better understanding of the 
screening in these materials |15j and a refinement of the Ohno potential are necessary. 
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Appendix A. Useful relations 

Appendix A.l. Peierls transformation 

The Peierls Hamiltonian is diagonalized via the transformation 

0,k,a = CtkCk,a + i/3fcCfc+7r,cr , b k ^ = f3 k C k ,a - iakCk+7T,a ■ (A.l) 

Its inverse reads 

Cfc.o- = a k a k ,a + Pkbk,a , quadc k +TT,a = -i(3ka>k,<r + ioffeafe.cr • (A. 2) 
Moreover, we have 

2t(5shi(fc) 2 2 2tcos(fc) 
Appendix A. 2. Help functions 

In order to express matrix second-order elements for the Hubbard interaction we 
introduce the abbreviations 

fi{fa,k 2 ) = a kl Pk 2 - Q!fe 2 /3fci = , 

f 2 (fa,k 2 ) = a kl ctk2 + f3k!l3k 2 = f 2 (k 2 , fa) , 

u 1 (k ll k 2 ,k 3l k i ) ee fi(k2,fa)h(ki,k 3 ) - f 2 (k 2 , fa) f 2 (ki, k 3 ) , 
u 2 (fa,k 2 , k 3 , fc 4 ) ee /i(fc 2 , fci)/ 2 (fc 4 , /c 3 ) + /i(fc 4 , k 3 )f 2 (k 2 ,ki) . (A. 4) 
The second-order matrix elements for the long-range Coulomb require the functions 
a , n 1 cos(fc) cos(fcr) . 

Vcos(A:) 2 + <* 2 sin ( fc ) 2 
„ . . 1 v— * <5 sin(fc) sin(fcr) „ , . , , 

i ^ ^cos(£;) 2 + J 2 sm(fc) 2 

-L/2<r<L/2,r|2=l ' ' ^ ' 

- cos(fc)sin(fcr)_B 5 (r)^j , (A.6) 
where (r|2 = 1) denotes all odd r. With the help functions 
C+(fa,k 2 )= ^— 2cos(r(fc 1 -fc 2 )) , 

C_(fci, fo) = y2 ~^TT~ 2cos ( r (ki - k 2 )) (A.7) 

we then introduce 

Vi(fa,k 2 ,k3,ki) = C+(fc 3 , fc 4 )/i(fc 2 , fci)/i(fc 4 , fc 3 ) 

- C_(fc 3 , fc 4 )/ 2 (fc 2 , fci)/ 2 (fc 4 , fc 3 ) ; (A.8) 
v 2 (fa,k 2 ,k 3 ,k 4 ) ee C+(fc 3 , fc 4 )/i(fc 4 , k 3 )f 2 (k 2 , fa) 

+ C^{k 3> fa)f 1 {k 2> fa)f 2 {k 4 ,k 3 ) . (A.9) 
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Appendix A. 3. Useful contractions 

We collect some contractions which frequently appear in the calculation of matrix 
elements to first and second order in the Coulomb interaction. The expression (0)o 
denotes the expectation value of the operator O in the ground state |FS) l|12[l of the 
Peierls insulator. The pair contractions are 

(nLa)o = \ , (A.10) 
(Wo = 5*r,o + 1 ~ ( ~ 1)T (Mr) + (-l) l B S (r)) , 

(ct a a k , a )o = -±=e- ikl (a k + i(-l)'&) , 
V L 

(*Jt tt ,)o =^ kl (p k +i{-l) l a k ). (A.ll) 
Important contractions with four Fermi operators are 

(niJt ,>fc 2 , CT )o = je i(kl - k2)l (h(k 2 , h) + i(-l)'/ 2 (fc 2 , h)) , 

+ i\ j e-^{-l) l {l-(-lY)a k f3 k , 

_i^ e ifcr(_l)I(l_(_l)r) afc/3fc . 



and 



S fcl ,, ( ft,,, - I) b+) = \e-'^-^ 1 h(k 2 M) 



2 J *Wn L 

+ i\e- i ^-W{-l) l h{k2,k l ) ,(A.12) 



+ 4e i(fel - fc2)i (-l) i /i(fc 2 ,fc 1 ). (A.13) 

Appendix B. Second-order matrix elements for optical excitations 

Appendix B.l. Second order in the Hubbard interaction 

As intermediate states we may have \m) = |11) and |m) = [21) . Due to spin symmetry, 
the contribution from |m) = 1 12) equals that from |m) = 1 21) . We find 



M 2 



St) ( e >P^P2) = (5a) I] p ,, : 

\ 2L ' klMMM e-E J= iE(k 3 ) 

[ S PiM ~ SpiM ± [^pi,fe 4 ~ ^Pi.fcsl] x {Pi -* P2} 

<5fc 1 -fc 2 +fe3-fe4,o["2(fcl, &2, &3, 

+ 4 1 -fc 2 + fe 3 -fc4:±7r[wi(fcl,fe,^3,fc4)] 2 , (B.l) 
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and 



M. 



U 2 \21). 



2,(a,t) (e; Pl ,p 2 )= [ - 



V 



E 



e-ELi^) 



£fcl-fc 2 +fc3-fc4+fc5-fc6!0 

|<Jpi,fci5Ai,fc 2 Wl(A;3, &4, &5, Afe) 

+ 5pi,fc 3 5fe 3l fc4«i(Axi fe, fesj fee) 

- ^ 1 ,fc 1 5fc 1 ,fc 4 ui(fc 3 , fc 2 , fcs, fce) 

"I" ^fci— fc 2 +fe3 — ^4 + ^5— fe6,±"' 

{flpi,fci<Sjfei,fc a u 2(fc3, &4> &5, M 

~ S Pu k 1 ^k u k 4 u 2 (k 3 ,k 2 , k 5 , fc 6 ) 

- ^i,fe 3 5fe 3 ,fc 2 u 2(fcl, &4i fc.s, fee)} x I 



P2 



} 



(B.2) 
Pi -> P2} 

The six-fold sum can be reduced to a three-fold sum. Nevertheless, expressions 
like this show that higher-order terms in the perturbation theory cannot be handled 
numerically because they involve five and more particle-hole excitations. 

Appendix B.2. Second order in the long-range interaction 

The particle-hole excitation in |(s, t)) can be destroyed. Therefore, the state |FS) with 
no particle-hole excitations also contributes, 

KlS&Pi'P^ = ± l) 2 «o(piMp 2 ))/(2e) . (B.3) 

Next, the states |11) contribute 



1 fv 



2 V L 



E 



o~k 1 -k 2 +k 3 -k 4 ,o 

{hi,kv$k 3 ,k4. L ( s piM v o(h) ± <5pi,fe 3 «o(^i)) 

+ (fipi ,k 3 - $P1 ,fei ) v 2 (kl, k 2 , &3, &l) 
± (<5pi.fe 4 _ Sp u k 3 )V2(k3, fc 4 , ki,k 2 )} 

x {pi -^2} 

+ Sk 1 -k 2 +k 3 -k i ,±ir 

{(^Pi,fe 2 ~ <^>i,fci) v i(h, k 2 , fc 3 , ki) 
± (<5pi,fc 4 ~ S pit k 3 )v 1 (k 3 , fc 4 , ki, k 2 )\ 



x {pi -> P2} 



(B.4) 
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E 



i 



ki<k3,k2<k± 
£fci-fc 2 +fc3-fe4,0 



x L [5 Pukl v (k 3 ) + 8 pi . k3 v a {ki)] 2 
(SpiM - S pukl ) v 2 (ki,k 2 , k 3 ,k A ) 
+ (t> P i,fc 4 - 5 P im) v 2(fa, fc 4 , k\, fa) 

- (^pi.fea ~ &vite) v 2(fa, fa, fa, fa) 

- (<5 P i,fc 4 - $pi,ki) v i(fa, fa, fa, fc 4 ) j 

+ 4 

{(<5pi,fc2 _ ^pi.fci) &2, fe, fa) 

+ (<5 Pl ,fc 4 - 4i,fe 3 ) wi(fc 3 , fa, fa, k 2 ) 
~ {&piM - &piM) v i(fa, fa, fa, fa) 

- (S Pl ,ki ~ S Plt k 1 )v 1 (k3, fa, fa, fc 4 ) j 



(B.5) 



Excitations with three particles and holes equally contribute. We have two equal 
contributions from 1 30) and 1 03) , 

2 



2M, 



V 2 \30) , 



v 



1 



k 1 <k 3 <k 5: k 2 <ki<k e e 12j=X E (fa/ 



^fel-fe2+fc 3 -fc4+fc5-fe6,±(/-l)w 

/=1,2 

|4 5 ,fc 6 4i,fe5 (vf (fa, fa, fa, fa) - (fa <-► fc 3 )) 
+ h 3 M S PiM ("/(fa, fa, fa, fa) - (fa *-> fa) 
+ S kl M S Pi,ki (v/(fa,fa, fa, fa) - [fa <-> fc 5 ) 

- hukAiM ("/(fa, fa, k 5 , fee) - (fa <-> fa) 

- hi,ke S Puki ( v f(fa,fa, fa, fa) - (fa <->■ fa) 

- $k 3 M S PiM (Vf (fa, fa, fa, fa) - (fa <-> fa) 

- h 3 ,kedpi,k 3 (vf(fa,fa, fa, fa) - (fa *-> fc 5 ) 

- hsM^PiM (vf(fa,fa,fa,fa) - (fa <-> fa) 

- 4 5 i 4 <5 Pl ,ic 5 (vf(fa,k 2 ,k 3 ,k 6 ) - (fci «-> fc 3 ) 

x{pi^ P2 }. (B.6) 
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Moreover, we find from the states 1 21) and |12), 

2 



2M 



y 2 |2i) 

2,(s,t) 



(e;Pi,P2) 



E 



fel </C3 ,&2 <&4 ,^6 
^fcl-fc 2 +fc 3 -fc4+fe5-fe6,0 

|<J/!3,fc4^Pl,fe3 ?; l(^lj ^2, fe) fee) 

+ ^fci.fca^Pi.fei^lC^! &4, fc.5, fce) 

- ^fei,fe 4 ^i,fci t; i(fe, fe, fej fee) 

- 5fc 3l fc 2 ( ^)i,*3 u x(Ai) ^4) fc 5, fee) 
± Sk 5 ,k 6 S Pl ,k 5 vi(ki,k 2 , k 3 , fc 4 ) 

T <5fc 5 ,fc 6 V.fc5 U l( fc l> fc 4, ^3,^2)} 
x [pi ^p 2 j 

|<5fc3,fc 4 ^pi,fc3 U 2(fcl, fe, &5 5 ^6) 

+ S ku k 2 Sp 1 ,k 1 V2(k3,k 4: , k 5 , fc 6 ) 

- S ku k 4 Sp 1 ,k 1 V2(k3,k 2 , k 5 , k 6 ) 

- 5k 3 M^piM v 2{k\,k A , fc 5 , fe 6 ) 
± Sk 5 ,k 6 ^piM, v 2(h,k2, k 3 , fc 4 ) 

T ^fe 5 , fee ^1,^5^2(^1, &4, &3, fe)} 



X [Pi - > P2| 



(B.7) 



The (5-conditions reduce the six-fold summations to a numerically tractable problem. 
Appendix B.3. Second- order mixed interactions 

Only two terms add to the matrix elements in second order. The states |11) give 



M. 



\ / U. 



1 



fel> fe 2 ,fc 3 ,fe 4 e -E J '=i^) 

Sk 1 -k 2 +k 3 -k i ,Q 

{(^P2,fc 2 - 5p a ,fci) U2(ki, k 2 ,k 3 ,k i ) 

i (<5p 2 ,fc4 ^P2,fe3) ^2(^3,^4,^1,^2)} 

x |(^pi,fc2 _ S pit k 1 ) v 2 (kx, k 2 ,k 3 , fc 4 ) 
± (^pi,fe 4 - S Plt k 3 )v 2 (k3, ki, ki, fe)} 

+ {Pl ^P2} 
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{(^P2,fe 2 - 5 P2M )ui(k u k 2l k 3 ,k 4: ) 

± (^P2,Ai 4 - ^p 2: fc 3 ) Wl(fc3, &4, h, fo)} 

x |(^pi,fc 2 - S Pll k 1 )v 1 (ki, k 2 ,k 3 , k A ) 
± (^pi,fe 4 - <5 pi , fc3 )ui(fc 3 ,A:4,fci,fc2)| 



+ {Pi P2} 

The states |21) and 1 12} equally contribute and give 



(B. 



2M 2,(s,t) '(e;Pi,Pa) = ^-p- J 



y 1 

e-EL%) 



7 „ ^fel-fe2 + fe3-fe4+fc5-feB,±(/-l)w 

/=1,2 

|<5fc 3 ,fc4^2,fc3 U /( fc llfe, &5, M 

+ 5 kl ,k 2 Sp 2 , kl Uf(k 3 ,k 4: , h, fc 6 ) 

- Sk 1 ,k i S P2 ,k 1 Uf(k 3 ,k 2 , h, fc 6 ) 

- Sk 3 ,k2 S P2M u f(ki,k4,, k 5 , fc 6 )| 

X {4 3 ,fc 4 ^pi,fc 3 w /( fc l: fc 2, M 
+ £fci,fc 2 V.fci W / (^3) ^4, ^5, M 

- S kl ,kJpiM v f(k3, k 2 , k 5 ,k 6 ) 

- h 3 M^pi,k 3 v f(ki, k 4 , k 5 ,k 6 ) 
± (Sk 5 M S PT-M v f(ki,k 2 ,ks,ki) 

- Sk 5 M 6 S Plt kr Vf(ki, ki, k 3 , A; 2 ))| 



-{pi 



(B.9) 
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